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^3 ; Abstract 

We consider the semileptonic and nonleptonic D s decay modes to final states 
with 77 and rf . We use QCD sum rules to determine the D s — > r\ form factor 
/+, and a generalized factorization ansatz to compute nonleptonic decays. We 
propose a parameterization of possible OZI suppressed contributions produc- 
ing the t]' in the final state, compatible with current data; such a scheme can 



be further constrained improving the precision of the measurement of the D s 



decay rates, as expected by the ongoing experiments. 



1 Introduction 



The exclusive D s decays to final states containing rj and rf represent nearly 30% of the 
total decay rate of the D s meson. Therefore, D s could be a suitable system where to 
gather information on important aspects of the rj — rf phenomenology, namely the long- 
standing issue of the r\ — rf mixing. Moreover, D s can be used to further investigate some 
unsettled aspects of nonleptonic heavy meson decays, such as the anomalously large rf 
production observed in several heavy meson decay channels. Examples are B~ — > K~rf 
and D° — > K rf, the measured decay rates of which are substantially larger than what 
can be expected by naive theoretical calculations. 

The current experimental situation concerning the D s transitions to 77 and rf is sum- 
marized in Table [l] [Q, mainly using the results obtained by the CLEO Collaboration 
in the past few years [El. The experimental results are expected to be improved in the 



Table 1: Experimental rates and branching fractions of semileptonic and nonleptonic D s 
decays to final states containing 77 and rf . 



Decay 


mode 


T (10- 15 GeV) 


/3(10- 2 ) 


Dt~ 


> 7]£ + V 


34.5 ± 9.3 


2.6 ±0.7 


Dt~ 


■ rfe+v 


11.8 ±4.5 


0.89 ±0.34 


Dt~ 


' 777T + 


22.6 ±6.7 


1.7±0.5 


Dt~ 


'VP + 


143.3 ±41.2 


10.8 ±3.1 


Dt~ 


' rfir + 


51.8 ±13.3 


3.9 ±1.0 


Dt~ 


>rfp + 


134.0 ±37.3 


10.1 ±2.8 



near future, since the analysis of the D s system is an important item of the experimental 
program of the current hadron facilities, as well as of the e + e~ machines running at the 
T(4S) peak. 

The results in Table [I] have inspired several considerations. First, it has been proposed 
that information on the 7] — rf mixing could be obtained just considering the semileptonic 
decay modes. As a matter of fact, writing the hadronic matrix element governing the 
transition D+ —>■ 7]l + v in terms of form factors: 

< rf V ')\- S1 ,c\D s { P ) >= fl{q 2 ){p + p')» + fl{q 2 )q, (1) 



B(D+ -> rf£ + v) 



could be 



(q = p — p') and a similar expression for Dt — > rfl^v, the ratio , , r . 

B(D+ -> r}£ + v) 

used to access the 77 — 77' mixing angle through the ratios of the form factors f± (q 2 ) / f±(q 2 ) 



which are related to the 77 — 77' mixing scheme |3], || . In particular, information could be 
gathered on the mixing scheme in the flavour basis |J |6| , which consists in writing the 77 
and 77' states as combinations of \rj q >= —=\uu + dd > and |?7 S >= \ss >: 

|7/ > = cos0 g \f] q > -sin0 s \r] s > 

\i]' > = sin(p q \rj q > +cos<f) s \rj s > . (2) 

It has been shown []5| that in this scheme a single angle is essentially required, since 
|0s — 0g|/(</>s + 4>q) C 1, a result confirmed by a QCD sum rule calculation 0. Therefore, 
one can safely assume <ft s ~ <p q ~ 0; the most recent estimates of <ft give values close to 
40° 0, §J- In the flavour scheme, the semileptonic form factors relative to Df — > 7]l + v 
and D+ — > r(l + v satisfy the relation 

rnm\ =cot * ' (3) 

so that the possibility of a direct comparison with the results for obtained from the 
analyses of other channels involving 77 — 77' particles could be envisaged. The situation is 
particularly simple in the case of semileptonic Df decays to positrons or antimuons, where 
essentially only the form factors /+ (q 2 ) are involved. However, in order to pursue this 
program, one has to neglect possible contributions to the semileptonic decay amplitude 
from diagrams where rj and rj' are produced through gluon emission; we shall consider 
this problem below. 

As for nonleptonic decays, naive factorization, using the semileptonic D s — ► 77 and 
D s — > 77' form factors and the Wilson coefficients relevant for the transitions in Table [I], 
does not allow to predict all the branching fractions of Df —> r]^Tx + and Df — > r)^p + 
0. The same conclusion is obtained by analyzing the various decay channels in terms 
of transition amplitudes related by SU(3)p symmetry to analogous amplitudes for D 
decays |T0j, or accounting for some effects of the inelastic final state rescattering fllT]. 



In particular, the prediction for the rate of the decay mode D+ — > i]'p + is lower than 
the experimental measurement by more than a factor of two. This is disappointing: in 
Cabibbo favoured hadronic D s decays the final state contains a single isospin mode, thus 
ruling out possible interference effects due to the elastic final state interactions; moreover, 
the conservation of G-parity does not allow to include inelastic effects of intermediate 
states consisting of ordinary mesons in the Df — > r]'p + mode [1^|. Therefore, a different 



mechanism must be invoked to explain the enhanced 77' production. It has been suggested 



that the enhancement could be due to OZI suppressed diagrams with the rf produced by 
gluons and the cs pair annihilating to a charged W || . This mechanism would not affect 
substantially the 77 production, since the coupling of the gluons to 77 is estimated to be 
smaller than the coupling to 7/ |13|| . However, a mechanism of this type, violating the OZI 



rule, could also affect the semileptonic -D+ — > t/£ + z/ transition, spoiling the possibility of 
using the relation (|3]) to gather information on the angle <fi from the semileptonic decay 
rates. Moreover, these effects could be also present in other systems, namely in D decays, 
although in such cases the annihilation amplitudes are Cabibbo suppressed. 

The effects of the gluon production of the 77 and r/', although plausible, are not included 
in ordinary analyses since they are difficult to take into account in a quantitative way. 
Nevertheless, their investigation is of particular relevance, and we shall try to perform it 
in a phenomenological way. 

In this paper, we compute the form factor f+{q 2 ) relative to Df — > r/£ + z/, showing that 
the result is in agreement with the experimental measurement in Table |1|. On the other 
hand, assuming the standard value of the 77 — 7/ mixing angle together with the naive 
factorization, other results in Table [I] are not reproduced. Therefore, we adopt a general- 
ized factorization ansatz, fitting the relevant parameters from the experiment; moreover, 
we assume that the effect of the process producing the 77' through the annihilation of the 
cs pair numerically modifies the D s — > 7/ form factors. This enables us to investigate 
whether the experimental results can be reproduced by this assumption and how the re- 
duction of the experimental uncertainty can be used to test various consequences of our 
ansatz. 

2 QCD sum rule calculation of f+(q 2 ) 

Let us first compute the form factor f+(q 2 ) using a nonperturbative method, such as the 
QCD sum rule technique [0. We adopt the usual strategy of considering a three-point 
function: 

IV(pW) = i 2 f*xd^e^e*~(0\T[J2(x)Jr(0)J?'(y)]\0) 

= ILh^ + IL.^ , (4) 

with J5 = S775S the pseudoscalar quark density probing the strangeness content of the 
r/, J M = s7 M c the weak current inducing the c — > s transition, and J 5 s = ci^^s a quark 
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current having the D s quantum numbers. The momenta P and q are defined as P = p+p' 
and q = p — p', respectively. For the invariant function U + (p 2 , p' 2 , q 2 ) a double dispersion 
relation in the variables p 2 , p' 2 can be written down: 

n+ (p ,p ,q) = ^ds! ds 2 - ^ ^ (5) 

it 2 J J (si — p 2 ){s 2 — p 2 ) 

where possible subtraction terms have been omitted. The spectral function p(si,s 2 ,q 2 ) 

contains, for low values of s\ and s 2 , a double 5— function corresponding to the transition 

D s — > r\. Isolating such a contribution, and neglecting possible subtraction terms which 

we discuss later on, we can write: 

Afl(q 2 ) f Da M 2 1 r r p had ( Sl ,s 2 ,q 2 ) 



n + {p\p'\q i )= ( ^ 2 " J ™> ^r,:- +^ / ds, I ds 



■J2~\ 



(M|, - p 2 ){M 2 - p' 2 ) m s + m c n 2 J s o J s o ( Sl _p2)( S2 _ p , 

(6) 
In (H) we have assumed that the contribution of higher resonances and continuum of states 
starts from the effective thresholds s® and s 2 - The hadronic parameter A represents the 
matrix element: 

<0|J 5 W)>=-4 (7) 

while the projection of the J® s current on the D s state is given by the matrix element 

<0\J°°\D s (p)>^- /DsM ^ 



m s + m c 

The correlator (Q) can be computed in QCD for large Euclidean values of p 2 and p' 2 
by an Operator Product Expansion, expanding the T-product in (B) as a sum of a per- 
turbative contribution and non perturbative terms, proportional to vacuum expectation 
values of quark and gluon gauge invariant operators of increasing dimension, the vacuum 
condensates. In practice, only the first few condensates numerically contribute, the most 
important ones being the dimension 3 < ss > and dimension 5 < sgaGs > condensates. 
The QCD expression for 11+ reads: 

n+(p ,p ,q) = -z dsi ds 2: 



■K 2 J(m c +m s )2 J Ami {s X - p 2 )(s 2 - p' 2 ) 

+ n+° =3) < ss > +n + D=5) < sgaGs > +... . (9) 

Invoking quark-hadron duality, i.e. assuming that the hadronic and the perturbative QCD 
spectral densities give the same result when integrated above the thresholds s? and s°> 
we get the sum rule: 

MoMk nw) _ i Ua /wl (10) 



m s + m c [M 2 - p' 2 ) {M 2 Da - p 2 ) An 2 Jd (si - p 2 ) {s 2 - p' 2 ) 

+ n+ D=3) < ss > +n + D=5) < sgaGs > 



with 

J> ert ( Q c r?\ 



and 



n 



(D=3) 



+ 2s 2 



1 ffflc 

2 l tt' 

+ (m c + m s 



s 1 — m 2 c + m 2 + s 2 + 2m s (m c — m s ) (11) 

4\/(s! + s 2 -g 2 ) 2 -4sis 2 

si + s 2 - g 2 )(2si - m 2 + m 2 ) (si + s 2 - g 2 ) 2 (si - m 2 + m 2 + s 2 ) 



(si + s 2 -g 



2N2 



4sis 2 



(si + s 2 -g 



2N2 



4sis 2 



(12) 



m 2 m c m s m 2 s (m\ + m 2 — g 2 ) 



2 2 

mi mi 



■ ^) i . i -. ' _ tyt^iiy' 



It t 









n 



(D=5) 



1 yQm s (m c + m s ) 6m 2 (m c + m s 



24 






/3 



<Y"J ty ' 



(13) 



2(m c + m s ){2m 2 + 2m 2 — 2g 2 — m c m s ) Y2m c 12m 



rv»^ (Y*' £ 



iy£ iyf 






r> 



The variables r and r' are defined as r = p 2 - m 2 and r' — p' — m 2 s . The domain .D is 
bounded by the curves 

± [2m 2 (si + g 2 ) + A(m 2 - m 2 - q 2 )} 



2m, 2 



(14) 



± 



2m 2 



(2m 2 (si + g 2 ) + A(m 2 - m 2 - g 2 ))~ - 4m 2 m 2 (si - g 2 ) 2 



1/2 



where A = s x — m 2 + m 2 , and by the lines s 2 = s° and Si = s°. Eq. ([II]) can be improved 
by applying to both its sides a Borel transform, defined as follows: 



B[T{Q 2 )} = Urn, 



Q 



2\ra 



d 



HQ 2 



(15) 



♦oo,n-oo,it-=^ (n-l)P " e ' \dg 2 7 

where JF is a generic function of Q 2 . The application of such a procedure to the sum rule 
amounts to exploiting the relation: 

exp(-s/M 2 ^ 



B 



1 

[s + Q 



2\n 



(M 2 ) n (n-1)! 



(16) 



with M 2 a Borel parameter. The operation, applied independently to the variables — p 2 
and — p' 2 , improves the convergence of the series in the OPE in the r.h.s. of Eq. @ 
by factorials in n, and, for suitably chosen values of the Borel parameters, enhances the 
contribution of the low-lying states in the hadronic representation of the correlator n + . 



Moreover, since the Borel transform of a polynomial vanishes, the procedure allows to get 
rid of subtraction terms in the dispersion relations, which are polynomials in p 2 or p' 2 . 
Therefore, a final sum rule can be worked out, keeping only the contribution of the lowest 
dimensional condensates: 



m s + m c 



fl( q 2 )e- M k/ M ? e - M '/ M 



_ e -ml/Ml e -m 2 JM. 



2 < SS > 

2 



m r 



+ 



An 2 Jd 
(m c + m s 



m 2 {m 2 + m 2 



dsids2p+ rt {si, s 2 , q 2 )e 



J ie 



2M| 



m c m s 

2 2 

mi /mi 



m,. 



2M(Mi 



2 Vm 2 4 Mf- 



)]} 



+e 



-mil M\ -mil Ml < SgoGs > 



+ 



WsK + m s) m 2 (m c + m s ) 

Mi + Mf 

2(m c + m s ){2m 2 + 2m^ — 2q 2 — m c m s 



3M{Mi 



4m c Am s 

~mJ 



M 2 



(17) 



In the numerical analysis of 
0.8 < qq > with < qq >= 



J17D we used standard values of the condensates: < ss > 
'—0.23 GeV) 3 , and < sgaGs >= m% < ss > with m^ 
0.8 GeV 2 . The charm and strange quark masses were fixed to the values m c = 1.4 GeV [1 
and m s = UOMeV M, O. As for the ,D S decay constant, we used f Ds = 225 MeV p 



while for the parameter A we adopted the two-point QCD sum rule result A = 0.115 Ge\^ 2 
0. The obtained sum rule shows stability to the variation of the Borel parameter in the 
region 2.5 GeV 2 < M 2 < 3.5 GeV 2 and 1.6 GeV 2 < M 2 < 2.4 GeV 2 , with the thresholds 
Si and s° m the ranges s? = 5.9 — 6.1 GeV 2 and s° = 0.9 — 1.1 GeV 2 , respectively. The 
form factor f1(q 2 ), obtained in the range of momentum transfer < q 2 < 0.5 GeV 2 , is 
depicted in Fig. [I]; it can be fitted by a linear expression 

fl{q 2 ) = Aq 2 + B, (18) 

with A = 0.14 GeV~ 2 and B = 0.50±0.04. This expression is consistent, in the considered 



range of momentum transfer, with a polar form /+ (q 
pole M P ~ 1.9 Ge^. 






2 ) 
P 



with the mass of the 



In the following, we shall consider the form factor f+{q 2 ) as a theoretical input in a 
phenomenological analysis of D s transitions. 



t 



o + 



0.8 



0.6 



0.4 



0.2 




0.1 0.2 0.3 0.4 0.5 

q 2 (GeV 2 ) 



Figure 1: Form factor f + a ~" n {q 2 ) as obtained using QCD sum rules. The shaded region 
represents the theoretical uncertainty related to the variation of the input parameters. 



3 D s transitions to 77 and rf 



The form factor f+{q 2 ) computed above allows us to calculate the semileptonic Df — > 
7]l + v decay rate. It can also be used to analyze the nonleptonic modes D s — > ??7r + and 
D s — > r]p + if the factorization approximation is adopted. This amounts to consider the 
effective Hamiltonian 



H. 



eff 



Gp 

V2 



v^v ud d0*) + 



CM 



(sc) V - A (ud) v _A + h.c. 



(19) 



with (qiq2)v-A = 9i7m(1 — 7s) 92 and C\ and C<i Wilson coefficients, and factorize the V — A 
currents appearing in it. As for the modes with rj', we further need an input on the rj — rj' 
mixing, and we choose the angle <ft i n the flavour basis mixing scheme, with the value <ft — 
39° coming from the measurements of <p — > rf , ' , y |I7 |. In Table ^| we collect the resulting 
branching fractions obtained in the factorization approximation, using f n = 0.132 GeV, 
f p = 0.220 GeV, td 3 = 0.496 ps; the number of colours N c is fixed to N c = 3, and the 
values Ci(m c ) = 1.263 and C2(m c ) = —0.513 are chosen, corresponding to the results for 
the Wilson coefficients obtained at the leading order in renormalization group improved 
perturbation theory at fi — m c ~ 1.4 GeV, in corrispondence to a s (Mz) = 0.118. Using 
the form factor /+ in (|18D we obtain the branching fraction B(Df -^ f]£ + u) = (2.3 ±0.4) x 
10~ 2 in agreement with the experimental outcome reported in Table [I]; also the result 

8 



Table 2: Computed semileptonic and nonleptonic D s rates and branching fractions. Non- 
leptonic rates are obtained using naive factorization. The 77 — 77' mixing is described in 
the flavour basis, with mixing angle <fi = 39°. 



Decay 


mode 


T (1(T 15 GeV) 


/3(10~ 2 ) 


Dt- 


f]i + V 


30.3 ±4.8 


2.3 ±0.4 


Df- 


rj'e+v 


12.7 ±2.0 


1.0 ±0.2 


Dt~ 


777T + 


38.5 ±6.2 


2.9 ±0.5 


Dt~ 


VP + 


74.5 ±11.9 


5.6 ±0.9 


Dt~ 


f]'TT + 


33.2 ±5.3 


2.5 ±0.4 


Df~ 


v'p + 


30.7 ±4.9 


2.3 ±0.4 



B(D+ — ► Vjl^v) = (1.0 ± 0.2) x 10 -2 , obtained using Eq.([|), is within the experimental 
uncertainty quoted in Table [TJ. On the other hand, as one can infer by comparing the 
computed decay rates reported in Table [^ with the experimental measurements in Table 
p], the calculations of the nonleptonic modes do not fit all the experimental measurements, 
as already anticipated by previous analyses. 

In order to parameterize the deviation from the factorization approximation, as well 
as the possible role of the 77 and 7/ gluon production, we adopt a generalized factorization 
ansatz, consisting in substituting the combination of the Wilson coefficients a± — C\ + 
— with effective scale-independent parameters a\ in the factorized amplitudes. The 



coefficients a\ should be considered as non-universal, process-dependent parameters |JI8 
However, since in the decay modes D^ 



r]7T + , Df —> v' n+ i an d analogously D~ s 



VP n 



D+ — > r]'p + , the underlying process is the same, we assume only two process-dependent 
parameters to describe the deviation from naive factorization: a\* describing Df — > r]ir + 
and Df -^ 77V" 1 ", and a{J describing Df -^ r]p + and D^ 



r]'p + . 



As for the possible contribution of OZI suppressed diagrams producing 77 and 77', it is 

where G is the gluon field and G 



(') 



essentially related to the matrix elements (0| GG 
its dual. Several theoretical investigations suggest that (0\GG\rj) <^ (0| GG \rf) (l^, 0; 
therefore, we assume that such annihilation amplitudes mainly affect the D s transitions to 
rf . A simple parameteterizion consists in modifying the values of the parameters A, B in 



18), thus without affecting the shape of the form factor f+ . This seems rather reasonable, 
since the range of momentum transfer in D s — > 7/ transitions is rather narrow (q 2 ~ for 

7]'p and q 2 < (Md s — M^) 2 for D s — > rfiv), and a linear 



D s -> 77V, q 2 



M 2 for D s 



9 



expansion is a suitable representation of the form factors. Therefore, in the case of ?/, we 
phenomenologically represent the D s — > rj' form factor as 



f e + ff (q 2 ) = Aq 2 + B. 



(20) 



It is now possible to use the experimental data in Table [I] to fit all the parameters we 
have introduced, namely a\^ , a\J, A and B. From the decays Df — > r/7r + and Df — > r/p + 
we find that the values of a{J and a\J are bound in the ranges: 

a[ f J e [0.65, 1.04] af p e i 1 - 21 . L86 l > ( 21 ) 

to be compared with the value of a,\ obtained from the Wilson coefficients C\ and 

C (m ) 
C2: ai(m c ) = Ci(m c ) H Ar ~ 1.1. As for the decay mode D 



N r 



r]'7i + , it involves 



/+ (<? 2 ); however, only the value /+ (0) = -B is needed in the approximation M w 
allowing us to constrain B in the range 



\B\ E [0.70,1.45]. 
Moreover, considering the modes Df — > rj'l + v and D~£ 



0. 



(22) 






(X P 



r(D+ -► r/£ + z/) 



x„ 



/f(M) 



,e// 

1,P 



7/p + , we find that the relations 

(23) 



X. 



semilep 



{M Ds -M^,f 



dq 2 A(Mi.,JWjY) /+'(? 



3/2 r 



re/// 



(24) 



32vrM|, s 



A(M|, s ,M 2 ,M 2 ,) 3/2 and X, 



semilep 



GfV cs 

1927T 3 M|, s 



constrain the pa- 



rameters A and i? in selected regions of the (B, A) plane. These regions are delimited by 
two straight lines, from the datum on the nonleptonic Df — ► 7]'p + decay rate, and by two 
ellypses corresponding to the measurement of the semileptonic Df — > ri'£ + i> decay rate. 

Considering simultaneously the constraints, all the data in Table [IJ can be fitted if, 
in the (B, A) plane, overlap regions exist among the area delimited by the ellypses from 
Eq. fl24T) , the regions delimited by the straight lines from Eq. fl23"|) and the regions between 
the vertical lines from Eq. (p22|) . At the present level of accuracy of the experimental data 
in Table [l] such regions indeed exist. They are depicted in figure || and denoted as Di, 
D 2 , D 3 and D4. The regions Di and D 3 are defined, respectively, by the conditions: 



1.2 < B < -0.70 



(25) 



2.9 [0.44 - 0.585] GeV~ 2 < A < 1.3£ 



-1.62S + 1.224V1.57-S 



GeV- 



10 




B 



Figure 2: Bounds on the parameters (B,A) in (p0|) . The ellypses represent the curves 
obtained from Eq. (p4|) ; the dashed lines stem from Eq. (|2"3[) ; the two pairs of continuous 
vertical lines represent the bound ( p^) . The shaded areas and the dots indicate the regions 
of the parameter space satisfying all the constraints. 



and 



-1.38 



1.625 + 1.224\/1.57- B 2 



0.70 < 5 < 1.20 (26) 

GeV" 2 < A < -2.9 [0.44 + 0.585] GeV" 2 . 



On the other hand, the regions D^ and D4 are defined, respectively, by the conditions: 



2.9 



0.44 - 0.585 



0.70 < 5 < 0.755 (27) 

GeV~ 2 < A < 1.38 [-1.625 + 1.224^1.57 - B 2 } GeV~ 2 



and 



- 0.755 < B < -0.70 {21 

-1.38 [1.625 + 1.224^1.57 - B 2 ] GeV~ 2 < A < -2.9 [0.44 + 0.585] Gel/" 2 . 



Although it is expected and rather plausible, the existence of such overlap regions was 
not guaranteed a priori; it shows that we have chosen a sensible scheme to parameterize 
the decays in Table [l]. More important, we expect that an improvement in the accuracy of 
the experimental data on the D s decay rates would sensibly reduce the size of such overlap 
regions, and presumably, exclude some of them. Noticeably, already at the present level 
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of accuracy some interesting observations can be drawn. Let us consider, for example, 
the parameters in the regions D3 and D^- In both the cases the experimental branching 
fraction of the semileptonic decay mode D s — » rflv is reproduced. However, a prime 
difference is that in the region D 3 the parameters A and B are opposite in sign, while in 
the region D2 they have the same sign. This implies that the relation between the D s — > rj' 
and D s — > rj form factors in (|3|) cannot be satisfied by the parameters in the region D%. 
The same conclusion holds for the region Di. The opposite signs between A and B, as it 
happens in the regions D\ and D3, have an observable consequence in the spectrum of the 
semileptonic decay Df — ► T]'£ + v. in this case, a zero in the dT(D s — > r]'£u)/dq 2 distribution 
should be observed, as depicted in Figure |3|. On the other hand, in the case of parameters 
in the region D2 (and D4) a smooth decrease in the spectrum should be observed as in 
D — > Klv. At the present level of accuracy of the experimental measurements, no choice 




Figure 3: Semileptonic spectra of D s — > rflv. The green curve corrisponds to the param- 
eters (B,A) = (0.9,-2.9) G £>i, the blue one to (B,A) = (0.72,0.08) G D 2 . 

can be done between the two shapes of the semileptonic distribution. 

We can reasonably expect that improved data would restrict the allowed regions in 
the (B, A) plane. It could happen that they do not intersect any more, or that intersec- 
tion regions could be found with restricted extension, allowing a better determination of 
the effective parameters introduced in our analysis. The calculation of such parameters 
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remains a challenging task, and we do not attempt it in the present paper. However, it is 
worth outlining the theoretical framework in which the calculation could be carried out. 

Concerning the effective coefficients alj and alj, which take into account the de- 
viation from the naive factorization in the corresponding decay modes, their theoretical 
calculation would consist in a precise determination of nonfactorizable contributions. A 
step in this direction has been recently performed in the case of some two-body nonlep- 
tonic B decays, where the meson picking up the B spectator quark is light, exploiting the 
large value of the beauty quark mass |nj. In this case, it has been observed that non- 
factorizable contributions are of order a s or 1/to&, and a QCD factorization formula has 
been written for the nonleptonic matrix elements in terms of meson light-cone distribution 
amplitudes. A possible extension of such a procedure to charm requires the development 
of a reliable method for computing at least the first (process dependent) 1/rriQ correction. 
A different approach would consist in considering the corrections to the large N c limit, 



where factorization becomes exact p0 |: also in this case, however, next-to-leading 1/N C 
terms are generally sizeable, and one needs their actual calculation. Therefore, it seems 
worth attempting to gain information on the effective coefficients from phenomenological 



analyses, as done, for example, in [18 



As for the OZI suppressed diagrams producing the rj' through its coupling to the 
gluons, together with a weak annihilation of D s , a perturbative calculation could be 
carried out in QCD, in analogy with the calculation of the rf production in quarkonium 
decays |H], 22[ • The difference, in the present case, is that one has to account also for the 



gluon emission from a light (strange) quark, and one cannot exploit the fact that all the 
quarks involved are heavy, which justifies the application of perturbative QCD methods. 
The calculation, for small values of q 2 , produces an amplitude for D s — > rj'iv of the same 
form as provided by a linear q 2 representation of the D s — rj' form factor. An important 
ingredient in this perturbative calculation is the actual value of the two-gluon-77' matrix 
element < g{ki)g(k2)\rj'(p) > describing the vertex rj'gg for off-shell gluons. Such a matrix 
element is parameterized by a form factor F{k\ , A;|) whose value at k\ = k\ = is fixed 
by the QCD anomaly; as for the momentum dependence, various parameterizations have 
been proposed in the literature, thus providing different values for the effective parameters 
A and B introduced in our analysis, which in turn could correspond to various solutions 
for the spectrum shown in fig.^. One might notice some analogies with the analyses 
which explain the observed enhancement of the rj' production in B decays through the 
mechanism of gluon fusion [2v| . 
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All such considerations taken into account, we believe that our proposed scheme, where 
additional contributions are reabsorbed in the parametrization of the D s — ► rj' form factor 
and in a\J, a\J is useful from the phenomenological point of view, as a starting point 
for the investigation of the underlying dynamics, and could be extended to other cases. 

Before concluding, we want to mention a check of consistency. If we consider the decay 
mode Df — ► K°K + , which can be related to Df — ► i]7i + through SU(3)f symmetry, 
and describe the D s — > K° form factor by f+{q 2 ), together with fx = 0.160 GeV, we 
can estimate the effective parameter a\ K . The experimental measurement B(D+ — > 
#°iir+) = ( 3>6 ± L1 ) x 10 -2 produces o^ G [0.72, 1.14], i.e. the effective parameter a\ f K 
displays a significant overlap with the range determined for a\ I . In different words, from 
our analysis and assuming SU(3)p, we would be able to predict rather accurately the 
experimental datum for B(D S — »■ K°K + ). 

4 Conclusions 

We have presented a phenomenological analysis of the D s decays to final states containing 
7] and i]'. Since the theoretical investigations based on SU(3)f symmetry, FSI effects and 
standard rj — rj' mixing failed in simultaneously reproducing the observed branching ratios 
for all these decays, we have considered a possible role of annihilation diagrams, in which 
the rj' is produced through its coupling to gluons. We have proposed a parametrization 
of those effects in the D s — ► rj' form factor. As for D s — > r], we used a theoretical calcu- 
lation of the form factor f+(q 2 ) which corresponds to a branching fraction for the decay 
D s — > rjiu in agreement with data. A fit to all the available experimental results, adopt- 
ing a generalized factorization scheme for nonleptonic decays, is possible; it constrains 
the parameters in restricted regions that can be discriminated by making dedicated ob- 
servations, for example looking at the semileptonic spectrum of the D s — ► rj' transitions. 
An improvement in the precision of the experimental data on D s decays could support 
this scheme and be helpful in understanding the dynamics of the rj and rj' production in 
heavy meson decays. 
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